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Free vibration characteristics of a t hin -walled 
open cross section beam with unconstrained damping layer is 
studied in the present work. A typical example of a top hat 
section with viscoelastic layer at the flanges is considered. 
Transverse vibrations both in and out of the plane (coupled 
bending and twisting modes) of symmetry are investigated. 
Results are obtained for natural frequencies and modal loss 
factors in both these modes of vibration. Beams with simply 
supported and clamped end conditions are investigated. 
Extensive parametric investigations are carried out with an 
efficient algorithm based on two dimension Newton-Raphson 
technique . 



CHAPTER 1 


INTRODUCTION 


1 .1 Introduction 

Control of vibration is an important aspect in the 
design of resonant systems. This is necessary in order to 
avoid early fatigue failure and malfunctioning of such sys- 
tems. .Vibrations control is also of significance in view of 
human comfort especially in vehicular systems. 

The level of vibration in a system can be reduced 
in a number of ways. These include, among others, 

(i) reduction of the excitation, (ii) isolation of the source, 
(iii) introduction of vibration neutralizers and (iv) detuning 
and decoupling by proper design. For resonant systems, how- 
ever, increase in the damping is the most effective way of 
controlling excessive vibration. Several damping treatments 
have been devised and practiced towards this end. For 
example, damping tapes are applied to aircraft structures with 
a view to increase the fatigue life and reduce the interior, 
cabin noise. Mastic deadners serve the same purpose in 
automobiles « 

The fuselage of an aircraft is subjected to con- 
vected random pressure field from the jet noise and turbulent 
boundary layer. As a result, the fuselage vibrates and in 
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turn radiates noise into the passenger cabin. Cabin noise in 
the mid and high frequency regions can be adequately atte- 
nuated by skin damping tape or fiberglass insulation. How- 
ever, these methods are not very effective at low frequencies. 

A concept similar to that of the tuned vibration 
damper has been attempted [ 1 ] in the design of the aircraft 
skin construction to reduce the vibration levels at low 
frequencies. This concept has been referred to as ’intrinsic 
tuning'. For a typical skin-stringer combination, the fre- 
quency of peak response is very close to the natural 
frequency of the individual skin bay (fp), clamped along the 
stringers and simply supported along the frames. If stringer 
and panel dimensions are such that the natural frequency of 

the stringer (f ) concides with f , then the stringer can act 

s p 

as a tuned absorber. To increase the effectiveness of the 
intrinsically tuned stringer over a broad frequency range, it 
is highly desirable to introduce additional damping. Addi- 
tional damping in the stringer can be incorporated by various 
methods. The most common methods are 

i) to have unconstrained or constrained visco-elastic 
layers on the stringer flanges, 
ii) to use inherently high damping composites or 
alloys as the stringer material. 

Stringers are thin walled open sections beam type 
structures. No adequate theory is available for the vibra- 
tion analysis of such structures with layered damping 
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treatment. The purpose of the present thesis is to make a 
theoretical investigation of such problems. 

1 .2 Short Review of Previous Work 

Work on layered damping treatments was initiated 
by Oberst and Lienard [ 2 ] . They considered the problem of 
an unconstrained layer damping treatment. The overall damping 
capacity of the 'composite plate was expressed in terms of the 
loss factor, Kerwin [3] investigated the problem of a cons- 
trained visco-elastic layer where loss modulus in shear 
contributes predominantly towards overall damping. Ross, 
Kerwin and TJngar [ 4 ] derived a general theory for free har- 
monic wave motions of an infinitely long beam made of three 
layers . An expression for the overall loss factor as a 
function of frequency was derived. The theory was valid only 
for simply supported end conditions in the case of finite 
beams. Ditaranto [ 5 ] presented a general theory for free 
vibrations of sandwich beams with arbitrary boundary condi- 
tions. Complex natural frequencies ( 00 ) were obtained in 
the form = u (1 + in ) where n is the loss factor. 

Interestingly, it was found that the modal loss factor 
bears a relationship with w 2 which is independent of end 
conditions. Moreover, this relationship was found to be 
the same as that given by the previous work on infinite beam. 

p 

However, the values of <0 depend on the boundary conditions 



and as such the numerical values of q need to he determined 
separately for each boundary condition. 

Mead and Markus [6 ] considered the forced vibration 
of a three layer damped sandwich beam with arbitrary boundary 
conditions. They used the concept of forced damped modes for 
the response calculation. Markus et al t 73 also discussed 
numerical techniques for determination of modal loss factors 
and natural frequencies of such structures. The problem of 
damped sandwich plates with elastic constraints along one pair 
of edges was investigated by Mead [8] . He also investigated 
the free vibration of a damped sandwich plate on periodic 
supports [9]* 

All these works, mentioned above, consider beams of 
solid cross section for which the centroid and the shear 
centre coincide. For such structures the bending and tor- 
sional modes are decoupled. 

As already mentioned, stringers are beam type 
structures with open cross sections for which the shear centre 
and the centroid usually do not coincide. In such cases 
there will be coupling between the bending and torsional 
modes of oscillation. Gere and Lin [10] presented an approxi- 
mate method based on Rayleigh - Ritz formulation for deter- 
mination of the coupled natural frequencies. 
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The major objective of the present work is to 
provide a theoretical analysis for the free vibration charac- 
teristics of thin walled open section beams with unconstrained 
viscoelastic layer. As an example, the problem of a top hat 
section with damping layer on the flanges is considered. The 
base material is assumed to be nondissipative and the loss- 
factor of the viscoelastic material is assumed to be indepen- 
dent of frequency. Moreover the loss factor is assumed to be 
the same in the extensional and shear modes of deformation. 

The exact analysis of the free vibration characteristics of 
thin walled open sections is presented in Chapter 2. Numeri- 
cal results of natural frequencies are compared with those 
obtained by the approximate method of Gere and Lin tlOl. 

The analysis for the uncoupled bending vibration in 
the plane of symmetry is essentially that of the transverse 
vibration of beams with solid cross-section [ 2 to 9 ) • The 
natural frequencies and modal loss factors are obtained in 
close form in terms of geometric and material properties of 
the beam. 

The analysis for the coupled bending and torsional 
oscillations require, as a first step, the determination of 
the St, Venant's torsion constant and the warping rigidity 
of the composite beam. These quantities turn out to be 
complex quantities due to the presence of the dissipative 
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layer. The St. Venant torsion constant is determined by 
extending the membrane analogy of elastic sections to the 
composite section. The warping rigidity is derived in terms 
of the average warping function. 

The coupled bending and torsional modes are gover- 
ned by coupled eighth order differential equation. Assuming 
harmonic motion for the free vibration, the frequencies turn 
out to be complex. The real part is the natural frequency 
and the ratio of bhe imaginary part to the real part gives 
the modal loss factor. 

The standard procedure of determining the frequency 
results in a frequency determinant with complex elements. 

The zeroes of this determinant determine the natural frequen- 
cies and modal loss factors * The complex arithmetic adds a 
dimension of difficulty to the determination of roots of the 
frequency equation. This difficulty is overcome by using a 
two dimensional Newton - Raphson technique. Numerical 
results are presented for various geometric and material 

parameters involved. 



CHAPTER 2 


FREE VIBRATION OF OPEN SECTION THIN WALLED BEAMS 


2.1 Formulation of the Problem 


Figure 2.1 shows a typical beam cross section 
considered in this thesis; where 0 is the shear centre and 
C is the centroid. X axis is the longitudinal axis through 
0 forming a right handed coordinate system with X and Z 
axes. In the following analysis v and w represent 
deflection along Y and Z axes respectively. $ refers to 
the rotation of the cross section about X axis. The coor- 
dinate system xtC bas its origin at C. 

It can be readily seen that 


V C = V Q - ° Z * 

where c„ is the distance of the centroid from 
z 

the shear centre, the subscripts refer to the points under 
consideration. 


E I 

E I 


The eauation of motion can now be written as 1 10] 

if 

3 w 
TjT 


3 x 
if 

3 v 


a 2 ‘ w 

0 + p A Y 


= 0 


(2.2) 


3 X 


3 1 

a t : 


5 (v o ' c z f ) = 0 


+ P A 


(2.3) 
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E C 


3 !+ \p 


w 


3 x 


Tj7 


G G 


- — * - P A (v - 


s 3x 2 


3t £ 


C Z ^ °z 


+ pi -M 1 = o 
0 at 2 


(2 


where and are area moments of inertia of the 

cross section about rj and ? axes respectively. A, C g and 
C v are area of cross-section, Saint Venant's torsion constant 
and the warping constant. E and G are elastic and shear 
moduli of beam material. I is the polar moment of inertia 
of beam cross-section. 


Equation (2.2) is the classical Euler's equation 
for vibration of the beam for which the results are well 
known . 

Assuming harmonic motions for the free vibration 


v = V (x ) e d w ^ 
o o 

and ifi - \fi Q (x) e d w 

equations (2 .3) and (2.*f) can be transformed, using equation 
(2.1 ) to 


+ ( 


d 8 V ( 
G C 


E C„ 
w 


E 


G G s d& V o . _p_A 
cbf ' EI C 

2 d2 V o 
2 , — Q, . / 

h ' ’ dx 2 


P a < pi 

+ T ~C~ + E c 
w w 


s. y- — EL-A ) w 

A. *rp ~r / 


P I ( 
E C. 


x 


p A \ ,.3f 

T ' 


w 


E I 


to ’ V 


o 


1 + 

m2 d V o 

dx* 


0 


(2.5) 
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In terms of non dimensional quantities, equation 
(2.5) can be rewritten as 


a 8 V d 6 V d 4 v 

— ”5 + a -i — r° + n*ou E ° 

d£ 8 1 dC 6 n 2 d^ 


O' 


n 


d 2 V k 
“3 +fi n“4 V o 


( 2 . 6 ) 


The constants <Xj , etc, are given in Appendix 

A y and is the nondimen sional natural frequency. 

8 A_ x 

Assuming the solution V (x) = V Q e 

0 s=1 3 

for equation (2.6), the following characteristic equation is 
ob bained for x 


8 


+ a. 


+ a. 


n 


2 ,4- 


n 


+ a. 


fi 2 X 2 
n 


+ a 


4- 


n 


4 

'n 


= 0 


(2.7) 


Applying appropriate boundary conditions at the 
two ends of the beam the frequency determinant (order 8x8) 
in the form 

Det [ A] = 0 (2.8) 

is obtained. 

The elements of the matrix A for simply supported 
and fixed end beams are given in Appendix A . The coupled 
natural frequencies, ^n are obtained by using Newton - 
Raphson method. 


2.2 Results and Discussion 

Table 2,1 shows the first eight natural frequen- 
cies for the coupled vibration. The table also includes 
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the results obtained by the approximate method of Gere and 
Lin [103, which are found to be in close agreement. Table 
2.2 presents the natural frequencies in the vertical (uncou- 
pled) and in the coupled modes, nondimensionaliSed with 
respect to the same quantities. It is seen that some of the 
natural frequencies in the vertical mode are quite close to 
those in the coupled mode. This suggests the possibility of 
these modes being excited simultaneously. This behaviour has 
also been observed experimentally [11 J. 
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TABLE 2.1 

COUPLED NATURAL FREQUENCIES BY EXACT AND APPROXIMATE 

ANALYSIS 


h = 
t 

where ai and w are the natural 
n a 

frequencies by present analysis and 
approximate analysis respectively. 


2.9400 cm 
1 .77 80 cm 
2.9718 cm 
0.2032 cm 


w 


Q 


n 


n 


*£E I 


V 


_c_ 


PA 1 


Tf 



p A 1 


END CONDITIONS 


MODE 

| Simply Supported Ends 


| Fixed Ends 

| Natural Frequency by 
j Present Analysis 

I 

fi Natural Frequency)! Natural Frequency 
| by Approximate 0 by Present 

| Analysis SI t 1 0 ] | Analysis SI 

1 

9.40 

9.40 

10.13 

2 

17.92 

20.86 

26.81 

3 

20*89 

23.91 

46.23 


38.91 

47.99 

84.90 

5 

67.21 

7^.96 

126.98 

6 

81 .60 

84.92 

174.91 

7 

104.23 

109.67 

232.42 

8 

149.40 

194.12 

248.94 
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TABLE 2 .2 

NATURAL FREQUENCIES IN COUPLED AND VERTICAL MODES 
OF A SIMPLY SUPPORTED BEAM N ONDIMEN SI ONALI SED IN 
THE SAME MANNER 

±2 = 2 . 5400 cm 

1 1 = 1 .7780 cm 

h = 2.9718 cm 

where ft and ft are natural 
t = 0,2032 cm c v 

frequencies in coupled and ver- 
tical modes respectively. 


MODE 

5 Frequency in Coupled 
jj Mode, ft rt 
! 0 

0 Frequency in Vertical 

1 Mode, a 

5 v 

1 

8,61 


9.88 

2 

28.60 


39.51 

3 

33.27 


88.89 

4 

61 .46 


1 58.0*+ 

5 

107,26 


246.93 

6 

130,23 


355.58 

7 

166.37 


483.99 

8 

238 .t-t- 


632.16 


0 ), 


\J 




PA 1 


4 



U PA 1 



CHAPTER 3 


EFFECT OF UNCONSTRAINED DAMPING LAYER ON . 

GEOMETRIC AND ELASTIC PROPERTIES 

3.1 Shear Centre Location 

When the damping layer is attached to the flanges 
of the cross-section as shown in Fig. 3*1 the shear centre 
of the composite cross-section will he shifted. The new 
position of the shear centre is determined by the following 
analysis . 

Consider the bending of the section about the Z 
axis under the application of a shear load V in the Y direc- 
tion through the shear centre 0. AA and BB are two adjacent 
sections along the length of the beam, distance dx apart . 

1 and 2 refer to the base material and the viscoelastic 
layer respectively. Consider a section CO at a distance y 
from the neutral axis Z . The imbalanced axial force due 
to the differential bending moment is equilibrated by the 
elemental shear force dF given by 

z o + t 1 

dF - 1 1 f ^ °1B " °1A' ) dz "i 

2 o 

z +t^ +t 2 

. + 1 1 j ( ct^ - a^) dz (3.1 ) 

z o +t 1 

where a refers to the bending stresses, with the subscripts 
indicating the material and section. 




Fig, 3-1 Section of the composite beam 
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The bending moment at a section is resisted by both 
the base layer and the viscoelastic layer. Assuming the 
longitudinal strain at a section distant y from the neutral 
axis is the same for both, the portions of the bending moment 
resisted by each can be shown to be 


M, 


M, 


MEi I C 


E -| 1 5 + E p 

i d 2 
M E£ I ? 

E 1 1 C,' + E 2 2 


(3.2) 


(3.3) 


where M is the applied bending moment and , 1^ , 

"1 2 . 

the area moments of inertia about the neutral axis, and and 
E 2 are the elastic modulus of the base material and the complex 
modulus of elasticity of the viscoelastic material. 


From Eqns . (3.1 ) to (3.3) and using engineering 
bending theory it can be shown that the shear stresses in 
the two layers are 

V E 1 Q 1 

( Vi = 

V E* Q 2 (3,1 

(Txy)2 = h V 

"where and Q 2 are the static moments about the neutral 
axis of the areas above the section at which the shear 


stresses are acting 
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Considering the moment balance about a convenient 
point in the plane of the cross-section of all the shear 
forces the shear centre can be located. The final expression 
for the shear centre location is given in Appendix B . It 
should be noted that the shear centre position in this case 
is not only a function of the geometrical properties of the 
section but also depends on the material properties. Because 
of the complex modulus of the viscoelastic layer the location 
is given by a complex number. 

3,2 Torsional Rigidity of the Two Layer Beam 

In the case of thin walled open section beams, the 
applied torque is resisted by two different types of defor- 
mation. The first part corresponds to the torque balanced 
by shearing stresses due to pure torsion (St. Venant's torsion). 
The other part is the torque resisted by the warping stresses. 

The St. Venant's torsional rigidity corresponding 
to the torque resisted due to pure torsion for the composite 
section is derived by an extension of the membrane analogy 
method. 

Consider the equivalent membrane problem, namely 
that of the deflection of bhe membrane made of two materials 
of equivalent size (Fig. 3*2a) subjected to internal pressure. 
The membrane wiich is initially in the form of 0 = 0 will 
deform in the shape 0 = 0 (z, y) as shown in Fig^ 3*2 b* 
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Y 



Fig. 3-2a Equivalent Fig.3-2c Approximate shape of 
membrane deflected membrane 




Fig. 3-2 b Actual shape of deflected 
membrane 
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If 1^ >> and >> the cross section 
of the membrane can be considered as independent of y 
except near the end as shown in Fig. 3.2c. 0 1 and 0 2 are 

the deflected surfaces of the membrane of material 1 and 2 
respectively which also correspond to the Saint Venant ’ s 
stress functions. 


Let Q o be the uniform twist along the length, 
then using the membrane analogy 0^ and 0 2 are governed by 
the following equation, 




2 G 1 ®0 



where 



0 

o 


and 



are respective shear moduli. 


(3.?) 


^ T xy 


The shear stresses can be determined from the stress 

9 0o 


functions using the relations 

8 0 

and (x xy ) 2 


J 1 


8 z 


3 z 


(3 .6 ) 


0 ^ and 0p can be found by integrating Eq. (3.5’) 
and applying the following boundary conditions'. 


0 1 ( 0 ) = 0 2 ( 0 ) 

0 ^ (* t.j ) = 0>2 ( ^2 ^ = 


At z 

i.e. 


= 0. 


( Vi 




3 0 


1 


1 


1 


3 z 


z=0 


G, 


* 


3 0 2 

3 z 


z=0 


where 


y represents the shear strain. 


(3.7) 



20 


0 2 


The final forms of 0^ and 0 2 turn out to be 
= " ° 0 (z 2 - t 2 ) + U G 1 © Q (z + t.j ) 



Q 

o 


where U 


(z 2 - tf) + vG* 0 q 

°2 <f - G 1 A 


G 2 * a + G i i 



( 3 . 8 ) 


The torque resisted by the shear stresses in the 
flanges induced by pure torsion is given by 


T f = 2 j J 0 dz dy 

= 2 / 0 1 d A 1 + 2 / 0 2 d A 2 (3*9) 

A 1 f A 2 

where A^ f and A 2 are the areas of cross-section 
of the flanges for material 1 and 2 respectively. The 
torsional rigidity of the entire composite beam including 
the torque resisted by the rest of the cross-section is 
given in Appendix B 


3,3 Warping Rigidity of the Two Layer Beam 

The warping stress generated due to fibers non- 
uniform torsion, induce shearing stresses which results in a 
re sis tive torque . 

The warping displacement w^ can be written as [17] 


w 


d 




M- 

d x 


( 3 . 10 ) 
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where 


w g = Average warping function, 
w g = Warping function, 
p - Rotation of cross-section. 


Since $ varies along the length of the beam, 
adjacent cross sections will not be warped equally and there 
will be axial strain e in the longitudinal fibres of the 
beam, given by 
3w, 


'x 


3x 


( w 


w„ ) 


d 2 lp 

dx 2 


(3.11 ) 


Assuming no slipping at the interface, the axial 
strain should be the same in both elastic and visco-elastic 
parts. But the warping stresses produced due to the axial 
strain are different due to different elastic moduli of the 
two materials . 


Assuming there is no lateral pressure between the 
longitudinal fibres, one can write 



Figure 3.3a shows an element cut out from the 
flange of the elastic portion of the beam shown in Fig. 3.3b. 
Along the middle line, the stresses vary with distance s 
measured from the edge of the section and can be calculated 
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from an equation of static equilibrium of the element 

m n o p. Projecting all forces on to the x - axis, one can 
wri be , 

9 ( t t ) ^ 

jj-j ds dx + t 1 — 3 -~ ds dx = 0 (3-13) 


Integrating with respect to s and observing that 
is independent of s and that vanishes for s = 0, 

the expression for can be writcen by using equation (3-12) 


r i h = ' E i v5 J ( W S - W S } h is 


(3.14) 


3x J o 

The resistive torque is obtained by summation along 
the middle line of the section of the moments of the elemental 
shear force about the shear centre. Therefore, the expression 
for the warping rigidity for the elastic portion of the 
beam can be written, using equation (3.14) 


M 


1 





3 X' 


m s 

J < J 


( w - w ) ds } r ds 

O o I 

(3.15) 


where m is the perimeter of cross-section and r is 
the perpendicular moment arm about the shear centre . 


equation 

equation 


Integrating the last part of the integral of 

m _ 

(3-15) by parts and using J ( v - w ) t. ds = 0, 

O s S I 

( 3 . 15 ) can be written as 

3 3 m 2 

= - 3L y f ( w - w ) t i ds 

1 3X J o S 


M 


1 


~ E„ C 
1 w- 


3x3 


(3.16) 



2b 


where C, r is the warping constant for the elastic 

W 

portion of the team as given in Appendix B. 


Similarly, one can proceed to find the expression 


for Mg, the warping rigidity for the visco-elastic part, as 


14 , 


wnere 


* 

- E 0 C 
2 w 


33 




3 


(3.1 7) 


Wr 


2 ax- 

is the warping constant for visco-elastic 


portion of the team and is given in Appendix B. 


Table 3*1 presents some calculated values of the 
quantities discussed in the Chapter for a Typical Cross- 
Section . 
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TABLE 3.1 

SHEAR CENTRE LOCATION, TORSIONAL RIGIDITY AND 
WARPING RIGIDITY WITH VISCOELASTIC LAYER 


~ = 0.001 , $ = 1.0 

h = 3.50 cm, 

E 1 = 70 x 10^ Kgf/cm 2 


=3*0 cm, 1 2 = 2.60 cm, 

t-j = 0.20 cm* 



0.50 

0.71677 

-O.OOO76 

0.10165 

0.10443 

0.70857 

0.17797 

1 .0 

0.71599 

-0.00154 

0.101 93 

0.38755 

0.70876 

0.36402 

1 .50 

0.71520 

-0.00233 

0.10249 

0.93828 

0,70896 

0.55827 

2.0 

0.71^1+0 

-O.OO313 

0,10340 

1 .8453 

O.70916 

0.76085 



CHAPTER 4 


FREE VIBRATION OF OPEN SECTION THIN WALLED BEAM 
WITH UNCONSTRAINED DAMPING TREATMENT 
AT THE FLANGES 


4.1 Free Vibration in the Plane of Symmetry 
If. 1.1 Equation of Motion 


Figure If»1 shows the beam cross section with a 
visco-elastic layer on the flanges b^ - b^ is the neutral 
axis of the untreated beam section, b 2 - b^ is that of the 
visco-elastic layer and b ~ b is the composite neutral axis. 


For displacement in Z direction the strain at 
b^ - b^ and b^ - b^ can be expressed as [^3 


— D 


p 

9 w 


9 x 


H 


9 2 w 


(4.1 ) 


20 9x 2 


where x is measured along the length of the 
beam and other quantities are explained in Figure 4,1. 


The net extensional force on each layer can be 
obtained from the strains at their respective neutral axes 
as 

2 

F. = - E. A. D 3 -— 

1 1 1 9 x^ 



jaJLw 

9 X 2 


and 


E 2 A 2 H 20 


(4.2) 
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As the total extensional force on the composite 


section is zero, one gets 


F 1 + F 2 ~ ^ E 2 A 2 H 20 - E 1 A 1 D ' ) 2 

d 


E 2 A 2 (H 21 “ D) 


E 1 A 1 D 


Thus the shift of the composite neutral axis from 

that of the original beam cross section, D, is given by 
„ * 

Er, A 0 H 01 

D = = D R + i D t (say) ... Of. 3) 

E 2 A 2 + E.j A^ n 

with E 2 = E 2 (1 + i g ), where 3 is the loss factor of 
the visco-elastic material, one gets 

E 2 A 2 H 21 (E 1 A 1 + E 2 A 2 } + 8 H 21 E 2 A 2 /k k) 

D-n “ * " - ‘ ~o O 'O "~'T vm- •*■+■ / 

R (E 1 A 1 + E 2 k 2 r + g^ e 2 a 2 


BH^ E^ A^ E 2 Ag, 


(E 1 A 1 + E 2 A 2 ) 2 + 3 2 e| A 2 


Of. 5) 


The equation of motion can be derived by using the 


relation 


3 2 m 


“ (p 1 A 1 + 9 2 A 2^ 


Of .6 ) 


The expression for the bending moment M can 


be written as 


M = + F 1 (-D) + M 22 + F 2 H 20 
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where ^ 


E 1 h 


3 2 w 
1 1 a x 2 


and M 


22 
with I 


* 

“a 1 * 


p,2 

d w 


( l . l . 7 ) 


TJ 


11 


22 3 x 

and I„ as the second moments area 
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of the layers about their own neutral planes. Using equa- 
tions (4.2), (4.6) and (4.7) the final equation of motions 
is obtained as 


B (1 + i v ) 


3^ w 


3 x 


¥ 


P A 



(4.8) 


and 


where 

B = 


P A = pLj A^ + P 2 A 2 
E 1 {I r^ 1 +A 1 (D E- D I )} +E 2 { 

A 2 ( + D R “ D I ” 2 D R H 21 " 2 P D I 



-h 2 i 

(4.9) 


)> 


VB 


2A i E 1 d h d i +2A 2 b 2 d i ®R- H aP + 

5E a { 1,^ + A 2 (H 2 , + D 2 -D 2 - 2 D r H 21 )> 

(4.10) 


Now, equation (4.8) is simply that of a hyst^ritically 
damped beam. So natural frequencies and loss factors at 
different modes for various boundary conditions can simply 
be written as follows: 
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For simply supported ends 

2 __ n^" tt^ 1 " B . _ _ 

w n ~ Tf pA ? n - 1, 2 ? 3? 

and n = v 

For clamped clamped ends 

2 - 2 ; 11 B 

w n 3[4 “ PA 

and n = v 


(*+•11 ) 


CI+.T2) 


If Eg « E.| , as is usually the case, B and v 
can be approximated from Eqs„ 0.9) and (4.1 0) as 


B * E 1 I n + E 2 (I n + A 0 H 2 ) 
1 'll 2 '22 2 21 

r2 


and n = v * g 


I n + A 0 Hi., 
E^ ] 22 2 21 


E. 


+ V E 1 {(I n 22 +A 2 H 21 )} 


(4.13) 


(4,1 4) 


4.1 ,2 Numerical Results and Discussion 

It can easily be seen from Eq. (4.11) and (4.12) 
that the loss factor in this mode of vibration in the plane 
of symmetry, is independent of both the end conditions and 
the mode number. 

Figures '4.2 and 4.3 show typical values of the 
loss factor for various values of tg/t^, Eg/E^ and the 
loss factor of visco-elastic material 3 . The following 
dimensions are used for the cross section of the elastic 



Loss factor (r\ ) 
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layer, 1^ =2.00 cm, 1^ = 2.60 cm^ h = 3*5 cm and 

= 0.20 cm. As expected from the approximate expression 
(4*14), the composite loss factor n is seen to be almost 
proportional to E^/E^ and also proportional to 0 * The 
nature of variation of n with t^/t^ is also as predicted by 
Eq* (4.14). 

4.2 Free Vibration in Coupled Bending Torsion Mode 


4.2.1 Equation of Motion 

Consider the top hat section with the damping 

layer in its flanges as shown in Figure 4.4. 0 and C 

respectively are the shear centre and the mass centre of 

composite section. Let v c and v Q denote the displace 

monts of C and 0 along the Y axis respectively. It 

is obvious that 


where c is the distance of the mass centre from the 
z 

shear centre. Considering the equilibrium ox elastic and 
inertia forces in Y direction in one hand and the moment 
equilibrium about a longitudinal axis through the shear 
centre on the other, the following equations of motion 

are obtained* 


(E 1 h. 


+ E 2 I ? 


4 

S v 




2 

3 v 

IT + P 1 A 1 + P 2 A 2 ^ ' ° 


i 


(4.16) 



35 




3x 


¥ 


C. 


9 2 ij i 

ax 2 


— ( p-j + p 2 A^) 


3 2 v 

J c 

st 2 z 


+ ( p, + p 2 i Cg ) 
❖ 


3 2 ^ 

3t 2 


= 0 


(4.17) 


where and are Young's modulus of elasticity, 

and I 5 are area moments of inertia about t> axis, A-, 

2 

and Ap are areas of cross-section, and are mass 

densities, 1 and I are polar moments of inertia about 
C 1 °2 

a longitudinal axis. Subscripts 1, and 2 refer to the base 
material and viscoelastic layer. , C 2 are the complex 
torsional rigidity and warping rigidity respectively of the 
composite cross section as derived in Chapter 3. 

Assuming harmonic solution of the form 


v = V (x) e lwt 
o o 

^ = *l > 0 (x) e iut 

and using the relation (4.15) eliminating , the equa- 
tion of motion can be transformed into the eight order 
ordinary differential equation in nondimen si on al form as 


d 8 V 
d $ 


o 


d 6 V, 


+ a . 


d V 


+ sr 


a 


4 

d V 
_ J_c 

2 


s? 



d S 2 


+ 



01 4 



(4,18) 


where a 1? a 2 , c^ } a k and SI \ are given in Appendix C. 
The final form of equation (4.18) is exactly similar to 
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the equation (2.6) for the beam without the layer. But 

1 ? a 2 ? “35 a L|. an<i are complex quantities in this 

case due to the presence of viscoelastic layer. In order 

bo avoid confusion, let these complex natural frequencies be 
* p p where 

denoted as 0 ^ (1 + i n ) /the real part ^ is referred 

c 

to as the natural frequency and n is the modal loss factor. 


4-. 2 . 2 frequency and Loss Factor Determination 


Assuming a solution of the form 


8 

V n (x) = l 

0 s =1 


Q e 
s 


X ,x 
s 


and substituting into the differential equation the follow- 
ing characteristic polynomial equation for X i s obtained 


R k 

X + Oj X® + 


-via- p 

■X -t- ot f 2 X^+a, n =0 
3 c 4 c 


(4.19) 


Applying appropriate boundary conditions at the 
two ends of the beam the frequency equation is obtained in 
the form of 


Let [A] = 0 


(4.20) 


where the elements of matrix A are given in 
Appendix C. 

The solution of the determinental equation in 
this case is much mote complicated as compared to that 
without the viscoelastic layer. This is so because of 
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_ 3Aj^ 

jAi R = 9~x J ~ ;T~u : and so on). The subscript j 

x y a •> 

indicates that the concerned functions are ev al uated at the 
point (x, The required partial derivatives in the 

iterations are computed by a central difference scheme. 
Starting from the lower end of frequency scale the itera- 
tions are continued to obtain the natural frequencies of 
various orders successively. The convergence criterion for 
each natural frequency and loss factor are specified in 
terms of the percentage change in the values in consecutive 
iteration and the absolute value of the determinant. 

Hie above scheme converged in most of the cases 
within six or seven iterations. The initial starting point 
for the variables did not present any problem in almost all 
the cases. However the determinant values vary over a wide 
range, with frequency. This presented some numerical prob- 
lems of dealing with very large numbers . This problem is 
alleviated by suitably scaling the determinant value at 
the bogining of the iteration scheme for the determination 
of a particular modal frequency. 

4-, 2. 3 Numerical Results and Discussion 

Numerical results are obtained, using the analysis ' 
presented in 4.2.2, for the following dimensions and pro- 

t 

perties of the elastic layer: 

1 1 = 2.00 cm, 1 2 = 2,60 cm, ^ = 0,20 cm and 

h s 3.^0 cm, P 1 = 2.7x1 0" 6 >Cgf - sec 2 /cm^ 

B 1 = 70 x 10^ Kgf /cm 2 , G 1 = 28 x 10’ Kgf/cm 2 , 

1 = 62.00 cm. 
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Table 4.1 presents the first eight natural 
frequencies and nodal loss factors for a set of typical 
values of the parameters involved. It is seen, unlike in 
the case of vibration in the plane of symmetry, that the 
loss factor depends on the end conditions as v/ell as modal 
number. This is due to the shearing of the viscoelastic 
layer caused by the coupling of bending and twisting modes. 
This fact has also been observed in the case of beams with 
solid cross-sections having constrained damping layer. 

For both types of end conditions the loss factor 
in the first mode is seen to be considerably higher than 
those at higher modes. No definite trend in the variation of 
loss factor is observed with the modal numbers. At this 
stage, no attempt has been made to have physical justifica- 
tion for this behaviour. Hie mode shapes are very complica- 
ted to visualize any pattern of deformation for the damping 
layer. Hence it is suggested to calculate individually the 
loss factor in the required mode . 

Figures 4.5 and 4.6 show the variation of loss 
factor with thickness ratio (^) for typical ratios of 

and ^2 and g f or simply supported and clamped ends. 

E^ P-j 

It is seen that in the case of coupled modes also the nature 

*b 

of variation with -r— is very similar to that observed in 

1 

the vertical mode (shown in Fig. 4.2). The modal loss 
factor also appears to vary linearly with 3 , the loss 

factor of the viscoelastic material. As already observed 
in Table 4.1, the loss factor in the first mode always 



TABLE 4.1 


NATURAL FREQUENCIES AND MODAL LOSS FACTORS IN 
COUPLED VIBRATION 


t, 

c 

t. 


E 


1.50, 


Pp 

~~ = 0.20 1 1 = 2.00 cm, X 2 


2.60 cm 


sf = 0 - 001 ’ 


P = 1 .00 t 1 = 0.20 cm, h = 3-50 cm 


The nondimensional coupled natural frequency 

2 


fl 


n 


0) 


E 1 I « 1 +B 2 I C 2 


where w is the natural frequency 


4 


( A^ + P 2 A 2 ) 1 
the composite loss factor. 


TsiMPLY SUPPORTED ENDS [ 


FIXED ENDS 


MODE 

NO 

a T 

I c L_. 

n 

I ft 

j c 

*r 

..i _n 

1 

3.85778 

0.00350 

7 . 56462 

0.001 72 

2 

13.37329 

0.00175 

20.23136 

0.00129 

3 

16.86532 

0.001 24 

39.20238 

0.00108 

4 

29.11707 

0.00126 

65.61231 

0.001 03 

5 

51.14014 

0.001 06 

103.62929 

0.00098 

6 

66.50152 

0.001 01 

133.62726 

0.00090 

7 

79.^5019 

0.00097 

177.67156 

0.00088 

8 

114.04927 

0.00092 

202.95039 

0.00096 



burns out to "be higher than that in the second mode. Com- 
paring the values of loss factors from Figs. 4.5 and 4,2, 
(for die same values of the parameters involved), it can 
be said that the unconstrained layer is more effecting in 
damping the lower order coupled modes. 

From Figs. 4.7 and 4.8 it is seen again that 

Ep 

the modal loss factor varies linearly with fr- . This 
trend was also observed for the vibration in the plane of 
symmetry (shown in Fig. 4.3). 

Figures 4.9 and 4.10 show that the loss factor 

Po 

is virtually independent of — in its feasible range. The 

P 1 

obvious suggestion is that one should go for light visco- 
elastic materials so that the increase in weight is minimum 
with the damping layer. 

Figures 4.11 to 4.16 show the changes in the 
natural frequencies with various parameters. 

It can be noted that the natural frequency is 

e 2 

virtually independent of 3 . For reasonable values of jjr- 

( - .001 ), the frequency is lowered with higher thickness 

ratio as the mass is increased without almost any change in 

E 2 

bending rigidity. However, with - .1, The increase 

in mass can be offset by the increase in rigidity so that 
the natural frequency increases with increasing thickness 
ratio. Figures 4.15 and 4.16 suggest the obvious results 
of lowering natural frequencies with increasing thickness 

ratio, the effect being more pronounced with higher values 
^ 82 
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Fig.4*14 Variation of natural frequency with thickness ratio 
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CHAPTER 5 


CONCLUSIONS 


An analysis for free vibration characteristics of 
a tain walled open section beam with unconstrained visco- 
elas bic layer at the flanges is presented. Transverse 
vibrations both in and out of tne plane of symmetry are 
considered. In the later case, the bending and the twisting 
modes are coupled. The natural frequency and modal loss 
factors, for both type of vibrations, are reported for 
various values of the parameters involved. 

The nature of variation of the modal loss factor 
with different parameters are found to be similar in both 
cases. The important observations include the almost 
linear variation of the modal loss factor -with the ratio of 
elastic moduli and with the loss factor of the viscoelastic 
material. However, in the coupled modes, the loss factor 
even with the unconstrained layer, depends on the end 
conditions and the modal number. The treatment appears to 
be more effective in damping out the lower order coupled 
modes . 
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APPENDIX A 


DETAILS FOR BIIN WALLED OPEN SECTION 
A-1 Sec bion Properties 

Various dimensions of the cross section are shown 
in Figure A-1 where 0 is the shear centre and C is the 
centroid. 

(i) Area A = 2 1^ t^ + 1^ + 2 h ^ 

(ii) Centroid position, C = { (1^ t^ + 2 h 2 t^ 

+ 2 1 1 t 1 (2 h + t 1 ) >/ 2A 


(ill) Area moment of inertia about r\ axis, 

1 t 3 

I = + | h 3 t^ + | 1 1 t 3 + 2 1^ t 1 h (h + t-j ) 


-AC 


(iv) Area moment of inertia about r, axis 


h H i 

~~ i — + 1 


+ \ t 1 + “V"" (1 1 + + 6 11 t 1 

h *b-j p 
+ ~2 1 2 


(v) Saint Venant's torsion constant 

~ „ 2 1 ,l 3 , 2 } t 3 , 1 1 t 3 

C s " 3 h h + 3 1 3 2 


(vi) Location of shear centre 


*1 ^ , X 1 * 2 * 

___ ( 2 + 2 it 

2 t 1 h , 1 2 , 1 1 n 


“X"" 


X + “3 




(vii ) 


Q + £ 


(viii) Warping constant 


°w “ b 1 { I i + x 2 + I 3 + \ + - n 

nhere m = 2 + l + 2 h 

- , h 0 *1 _ 

¥ s = ~2 1-j Oh + e) g- — 

1 - 1 = (h + e ) 2 lj / 3 


_ 2 

m vr } 
s 


lg " T 2 { 12 (h + e ) 2 + h 2 l 2 

- 6 h lg (h + e) } 

1 3 ~ ^pg ^ 3 s + 4 e lg - 6 lg e a } 


where a = 1 2 h - 2 L (h + e) 

: l( . - “g { 3 h 2 + 3 h lg h + h 2 l 2 } 

where b - lg h - 2 1-j (h + e) - 2 e lg 

: r B -y-{ 3 a 2 + (h + e) 2 l 2 + 3 1 1 (h + e) a > 

where d = 1^ (h + e) + e lg ~ h lg 


A- 2 Won Dimensional Terms in Equation (2,6) 


(1 + p.j + s.| ) 


' i2 
§•] 1 


and a 



where 


G C 

rrj S. 

b 1 EC 

'VI 

T «2 

p ' - 

Pl " G w 

i i c T t 

CJ — >L V 2. 

***“ /“X A A 


1 is the length of the beam. 


Non dimensional natural frequency 


o 2 - “ 

n El 

__ -C^ — 

p A 1 

where w is the natural frequency 


A. 3 Appropriate Boundary Conditions for the Coupled Vibration 
For simply supported ends the boundary conditions 


at x = 0 


and 3c = 1 


and 


1|> ~ 0 
. ii rv 


- o 

o 

= 0 

0 


For fixed - fixed ends, the boundary conditions are 


at x 


r iv = 0 


and x = 


= 0 



A- 5 

A.*+ Elements of Matrix A in Equation (2.8) 


Elements of the matrix A is different for diffe- 
rent boundary conditions. For simply supported endss 



. 





1 .0 

1 .0 

1.0 

... i.o 1 


h 

e 

e * 2 

e * 3 

... e 8 


A 1 

A 2 

X 2 

A 2 

3 

... *§ 

[ A ] = 

, 2 Al 

A -j © 


A 2 e^ 

3 

- ^ 


A 1 *" 

A 1 

A Lh 

2 

A^ ♦ * • 

3 

X^ 


A -j 6 

\ 4- A 2 

A 2 e 

X^ e X 3 ••• 

3 



A 6 

A 1 

A 6 

A 2 

16 

A 3 

X 8 


,6 *1 

A,. e 

6 ^2 

Ao e 3 ... 

... X 8 e * 8 


Here, X| , A . .. Ag are eight roots of the characteristic 
polynomial equation (2.7). 



For Fixed ends beam: 


X^ e 


X^ e 


X~ e 


X g e 


h 0 


F e X2 
2 e 


X^ e * 3 


X 4 - a X S 
X 8 e 


r 5 « *1 


xs © 


3 \ 

4 e 


*3 e " 3 


3 . 8 



APPENDIX B 


GEOMETRIC AND ELASTIC PROPERTIES TOR COMPOSITE 
CROSS SECTION 


p.1 Sec tion Properties 


re B.1 shows the composite cross section 

whore tl - R, and n, - are tlle centroida:1 - axis of P aI “ ts 
1 and 2 respectively. C is the mass centre of the 
composite section, 

J?I 5 * A 2 P 2 \ 

q ~~ A 1 P-j + ^2 

where are given as A and 5 respectively in 

Appendix A. and p 2 are the mass densities of the 
material of 1 and 2 respectively. 

t 2 

A 2 ~ 2 1-1 (h +t ! + "2 ) 


~ ^ given in Appendix A. 


j t 2 + a, + id' 


where I ? and I ? are the area moments of inertia 
of part 1 and 2 about the axis of symmetry . 


13,2 Shear Centre Location for Composite Beam 

I 0 

^ — \/ _ I V -4 ** ’ ^ 


V 3 1 2 - V 1 (2 h + t 1 + t 2 ) - 2 h 




whore 


B-3 


V, 


E i. { +1 2 } t i + 1 1 i h 2 > 


(E„ I 


-i !>. + E 0 I ) 

? -i 2 £' 


V, 


and V 


1 


with E, 


* 


t 1 E 1 _ 2 1 p li 

(•S-j Ir + E* I ) 1 t ^ + ~3 } 

1 2 c 2 

^2 E 2 2 / X 2 X 1 

«. v’^v 1 

E 2 ^ + 1 5 &= loss factor of the 

viscoelastic material. 


L«3 Coruional Rigidity of the Composite Beam 


C -j s 1 * 1-j ^ ( | ~ y t 2 ) 

+ k 1 1 G* ( | t| - 1 p t|) 

+ 0.. C ^ h + ^ Xg ) 

Gp t; - g 1 tf 
whom y n — - 

G a h + G i h 

with Or* = tfg (1 +i e ) 

B.'f Warping Rigidity of the Composite Beam 

°2 a E 1 t 1 G w 1 + E 2 h C w 2 

where C = warping constant of cross section 1 
as C,_ in Appendix A . 

W 



where to 


( 



APPENDIX C 


DETAILS POD THE FREE 'VIBRATION OF COMPOSITE 
BEAU IN COUPLED MODE 


C.1 Section Properties 

The various section properties can he written 
for the composite cross section shown in Figure B-1 , as 

2 


I ~ I + A.,(5 rt -C-,) 

"i 1 


•n< "n L i ^ c -1 


.3 


TV 


li 

+ 


1 2 ' A- ( ?„ - Co) 


X = I + I 
°1 X n 1 I?1 
Ic 2 = ”2 + C 2 

where !„ , I c , \ ^ C 0 . C 2 etc - are « iven 

in previous appendices , 

0,2 Non Dimensional Terms in Equation (4.18) 


a. 


a. 


a 


a 


4 


where 


■ g 1* 

• ( q + P + s) 

g q i 2 

s q 

C„ 


C and. C 0 are the torsional and the warping 
1 2 

rigidity given in Appendix B . 



c-2 


cl 


( 


p a . n a •S-t I* + Eq I r 

:1 1 p 2 a 2 ) x( J c i 2 


E- X r + E* I r 

' i 2 C 2 


P t A 1 + P 2 A 2 
>!< 

where E 0 is the real part of E„ . 

r- d 


p .. A., + P 9 A P p E 1 I ^1 + E 2 \r> 

( — J- - - L. — i — - c? ) x ( r-arJ—r— ) 


C 


P 1 A 1 + p 2 A 2 


P 1 I c 1 + p 2 I Cp E 1 + E 2 I ^ ? 

3 “ ^ C 2 ^ x ( ? 1 “'Tj" + p 2 'X 2 “ ^ 

Non dimensional. natural frequency is given hy 


o 


n? (i + iv) = 

* 1 c 


c 


b 


( p ^ p 2 A 2 ^ 2 

where w is the natural frequency 

ft is the nondimensional natural frequency 
c 

and n is the composite loss factor. 

v « n . 

0.3 Elements of Matrix A in Equation Of .20) 


Thu f orm of matrix A is same as given in 
Appendix A.4. 'lore, ^ ^ ... are the eighth root of 

characteristic polynomial equation given in Equation (4.13) 
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